Astronomy & Astrophysics manuscript no. 0615vb 


February 2, 2008 


(DOI: will be inserted by hand later) 





A numerical model for the 7-ray emission of the 
microquasar LS 5039 

V. Bosch-Ramon, J. M. Paredes 
' Departament d'Astronomia i Meteorologia, Universitat de Barcelona, Av. 



o 

(N 



Diagonal 647, 08028 Barcelona, Spain 
e-mail; vbosch@am.ub.es; jmparedes@ub.edu 



the date of receipt and acceptance should be inserted later 

, Abstract. The possible association between the microquasar LS 5039 and the 

, EGRET source 3EG J1824— 1514 suggests that microquasars could also be sources 

, of high-energy 7-rays. In this paper, we explore, with a detailed numerical model. 



o 



if this system can produce the emission detected by EGRET (>100 MeV) through 



. inverse Compton (IC) scattering. Our numerical approach considers a population 

O ■ 

of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting 
with both the radiation and the magnetic fields, taking into account the Thomson 
Q ■ and Klein-Nishina regimes of interaction. The computed spectrum reproduces the 

■J-J ' observed spectral characteristics at very high energy. 

Key words. X-rays: binaries - stars: individual; LS 5039 ~ gamma-rays: observa- 
' tions - gamma-rays: theory 

w. 

1. Introduction 

Microquasars are a selected class of X-ray binaries that produce relativistic radio jets 
(Mirabel & Rodriguez II 9991 Fender l^004|l . The origin of the jets is related to the matter 
accreted by the compact object, a neutron star or a black hole, from the companion star. 
These systems behave as scaled-down versions of quasars and active galactic nuclei. The 
population of microquasars is still a very reduced one, with about sixteen known objects 
up to now (Ribo 120021). Several models have been developed to explore the high energy 
emission from the jets of microquasars. Two kinds of model can be found in the literature 
depending on whether hadronic or leptonic jet matter dominates the emission at such 
an energy range: the hadronic jet models (e.g. Romero et al. I2003|l . and the leptonic jet 
Send offprint requests to: V. Bosch-Ramon 
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models. Among leptonic jet models, there are IC jet emission models that can produce 
X-rays and 7-rays, based in some cases on the synchrotron self-Compton (SSC) process 
(i.e. Band & Grindlav 119861 Atoyan & Aharonian 1 1 999() . and in other cases on external 
sources for the IC seed photons (EC) (i.e. Kaufman Bernado et al. 121)021 Georganopoulos 
et al. 12002(1 . In addition, there are synchrotron jet emission models that can produce 
X-rays (i.e. Markoff et al. lMIUl . 

The microquasar nature of LS 5039 was clearly established when non-thermal radi- 
ation produced in a jet was detected (Paredes et al. I2000|l . Interest in this source has 
grown significantly because it has turned out to be a source of persistent radio emitting 
relativistic jets (Paredes et al. 12002(1 . and because of its proposed association with the 
unidentified source 3EG J1824-1514 (Paredes et al. lWMlll of the 3rd EGRET catalogue 
(Hartman et al. 11999(1 . Recently, CoUmar 1(2003(1 reported the detection of an uniden- 
tified 7-ray source by the COMPTEL experiment, and LS 5039 is one of the possible 
counterparts. We are interested in investigating if microquasars are able to generate the 
high-energy 7-ray emission detected by EGRET^, using a numerical model based on IC 
scattering, taking into account energy losses, EC and SSC interactions. In this paper, 
we focus our attention on the microquasar LS 5039. In Sect. [21 we present the physical 
scenario of our model, and the model itself is described in Sect. El We apply the model 
to the microquasar LS 5039 in Sect. 01 and the obtained results are discussed in Sect.O 

2. The physical scenario 

The scenario of our model consists on an X-ray binary system where the compact object 
has a disk coplanar with the orbital plane and a jet not necessarily perpendicular to 
it. It includes also a hot spherical region, called corona, surrounding the compact ob- 
ject. We assume that the leptons dominate the radiative processes that generate 7-rays. 
Relativistic electrons, already accelerated and flowing away into the jet, are exposed to 
photons from the star (assumed to be an isotropic photon field) as well as to photons from 
the synchrotron emission of the same population of electrons, since we suppose the pres- 
ence of a magnetic field. This magnetic field is tangled in such a way that the resulted 
synchrotron emission is isotropic in the jet's reference frame. Although some features 
like disk-jet interaction (e.g. for blazars, Dermer et al. 11992(1 . reacceleration of the jet's 
particles (e.g. Markoff et al.^UU^, corona-jet interaction and pair creation-annihilation 
phenomena (e.g. both for Cygnus X-1, Romero et al. 12002(1 could be significant at some 
level, their study is beyond the scope of this paper and they will be treated in upcoming 
works. The 7-ray emitting region, the 7-jet, is assumed to be closer to the compact object 
than the observed radio jets. 



^ http:/ /cossc.gsfc. nasa.gov/egret 
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3. Description of the model 

We have developed a numerical IC model for the 7-ray emission from relativistic electrons 
in a cylindrical jet (the 7-jet). This 7-jet is assumed to be short enough to allow the 
magnetic field (B^) to be considered constant. We have included the interaction between 
the relativistic electrons and both the magnetic and the radiation fields. The energy 
losses of the leptonic plasma in the 7-jet are due to its own synchrotron emission, the 
SSC scattering and the EC scattering of the stellar photons by the relativistic electrons. 
Due to the strong losses, the electron energy distribution density along the 7-jet model 
varies significantly, and in this sense the 7-jet is non-homogeneous. Thus, the 7-jet is 
studied by splitting it into cylindrical transverse cuts or slices. The size of the slices 
has to be suitable to get almost homogeneous physical conditions within them (energy 
densities for the radiation and the electrons). 

The code used for our 7-jet model is run in two steps. The first computes, for each 
slice, both the total radiation energy densities and the radiation energy densities per 
frequency unit from both the star and the synchrotron process. Then, the calculated 
radiation energy densities are introduced jointly with the electron distribution function 
(see Sect. 3.1) in the second step of the program to calculate the IC emission at each 
slice. Because such slices evolve due to the changes of both the electron energy and the 
radiation energy densities, the resulting IC radiation will be different depending on the 
age of the slice (i.e. depending on the distance to the compact object). Once we have 
computed the IC emission from all the slices, we can sum the contribution of all of them 
to get the whole emission. Since our model is applied to a wide range of energies for both 
the electrons and the seed photons, we have taken into account both the low and the high 
energy regimes (i.e. the Thomson and Klein- Nishina regimes). We have used the cross 
section of Blumenthal & Gould H197(J|I . The free parameters of the model are just three: 
B^, the maximum electron Lorentz factor (j^^^) and the leptonic kinetic luminosity or 
leptonic jet power (iko)- 

3.1. Energy distribution of the electrons 

The differential equation that describes the electron energy evolution in our model is: 



where ag and a^ are constants (of values 2.37 x 10""^ and 3.97 x 10^^ respectively, in 
cgs units), C/star is the companion star's radiation energy density, Ugyn is the synchrotron 
radiation energy density, computed using the local approximation (Ghisellini et al. ll985|l . 
and E is the electron energy. It is worth noting that we work now in the jet's reference 



dE_ 



(as B^ + Oc t/syn + Ac t4tar) 



(1) 
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frame. Solving this equation for the given conditions, after some algebra we get the 
electron Lorentz factor evolution function: 



7o 



7oO 



1 + 7oomoc2[Cs(t - to) + '1 h"- (arctan 

/cu c L V ^1 / ' .t' Sin -cj ^ t' sir 



arctan 



tp—t' cos OT ^ 



(2) 



sin zu -'-I 

where t is the time and to is the starting time in slice evolution. Electrons are injected at 
a distance zq to the compact object, which corresponds to the time to, and their evolution 
continues along the jet axis (to farther z). is the jet velocity, tu is the angle between 
the jet and the orbital plane, and 700 is the value for the electron Lorentz factor at to. 
i?orb is the orbital radius of the binary system, and t' is i?orb/i^jct- Cs is (flgS^ + acC/syn)j 
and Cc is acrjctistar/47rc, where Lstar is the bolometric luminosity of the companion star, 
Fjct is the Lorentz factor of the jet and c is the speed of light. We describe the energy 
distribution of the (injected) electrons, N , with a power law in the slice at to- The initial 
slice (like the following ones) is assumed to be spatially homogeneous and the particle 
motion direction follows an isotropic distribution (in the reference frame of the jet). So, 
at to'. 



iV(7co) = Q7c7. for 7ct" < 7cO < 7 



max 
cO 



(3) 



where Q is the normalization constant of the electrons, related to the jet power 
(Georganopoulos et al. l2()()2|l . and p is the power law index of the electron energy distri- 
bution function. 7"g'''^ is the maximum Lorentz factor of the electrons in the first slice, 
which will be treated later. 7^'" is the minimum electron Lorentz factor in the first slice, 
whose value is determined according to both the upper limit of the seed photon energy 
and the lower limit of the outgoing photon energy (see Sect. 3.3). Taking into account the 
conservation of the number of particles, N{'jg,t)dje — A^(7eo)c^7eOj the time evolution of 
Eq. Q is found: 



1 - 7omoC^ Cs(t - to) + 



vfett' sin TO 



arctan 



t — t' cos tn 
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arctan 
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(4) 



3.2. The radiation fields 



Assuming that the stellar photon energy density is isotropic, C/gtar in the reference frame 
of the jet is given by rjetistar/47rcd^, where d is the distance between the companion 
star and a certain slice of the 7-jet. Assuming that the star emits like a black-body and 
taking into account its spectral type (see Sect. 4.1), we can obtain the radiation energy 
density per frequency unit (C/^ star)- The numerical calculation of C/gyn is difficult when 
this seed photon field becomes the dominant one. The problem is to know previously 
the evolution of the Usyn within a slice along the 7-jet. To know Usyn, Af(7o,t) has to 
be determined, but iV(7c,t) depends on the value of J7syn in the previous slice, as can 
be seen in Eq. Q. We solve this in the following way: Usyn and the radiation energy 
density per frequency unit {Ui, syn) are found initially starting with Eq. Q, determining 



Bosch-Ramon & Paredes: 7-ray emission in microquasars 5 

their own values for the next step and so on. The synchrotron emission is numerically 
integrated following the standard formulae from Pacholczyk H1970|l . The two integrals 
for computing numerically the emission and absorption coefRcients are: 

j^{z) = C3B-y N{jc,z)F{v/vc)dja (5) 
kAz) = -7; / 3— 2 7c(^) F{v/vc)djc (6) 



where C3=1.87 x 10"^'^ (in cgs units) and F{v/vc) ~ x kt^/3{z)dz, in which K^/3{z) is 
the Bessel function of the second kind, v and Vc are the frequency and the characteristic 
frequency respectively. After some algebra, C/^ syn is obtained. 

3.3. IC interaction 

The code developed for reproducing the radiation-matter IC interaction operates over 
files created by the first step. The seed photon energies (cq) should be such that the 
scattered photons had energies inside the chosen energy band following the relationship: 
46o7g 

'"^^^ 1 + 4607c ^' 
In the second step, both the radiation and the matter interact throughout the cross 
section of the IC process (Blumenthal & Gould [T1)7(J|I : 

4eo7c 

where ctt is the Thomson cross section constant and f{x) is: 



(T(a;,eo,7e) = 7— ^ /(x) (8) 



2a; In a; + a; + 1 - 2a;^ 



(4eo7ca;) 



2 



2(1 + 4eo7ea;) 
in which x is: 



P(l/472, 1, x) (9) 



^ = 1 — 27r — (1"^) 

and 

P(l/47^^ 1, x) = 1, for 1/472 < a; <1, and otherwise. (11) 

P(l/47^, 1, x) restricts the cross section to a physical case, where e cannot be lower 
than eq or higher than emax- 

The resulting IC spectral energy distribution in the observer's reference frame is: 

eL, = Vsii,,{z) / " ^ V ^° " Utot.iieo, z)N{j,, z) a [x, eo,lc)djdeo (12) 

Here, (5 is the jet Doppler factor, viz. S — [Tj^t (1 ~ l3coa9)]^^ and /3 = fjct/c VgUcc is a 
slice volume, which varies depending on the accuracy needed by each slice, 9 is the angle 
between the jet and the observer line of sight, and J7totai(eoi z) — C/syn(eo, ^) + C/star(eoj z), 
which is the total radiation energy density per seed photon energy unit and slice, taking 
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into account both the synchrotron and the star radiation field components. The integral 
has been calculated numerically taking suitable integration steps from the behavior of 
the involved functions. We note that this is the optically thin case because, in the 7-jet 
conditions and at the EGRET energy range, the IC absorption coefficient is negligible. 
It is worth mentioning that the electron energy evolves following Eq. |(5J. Therefore, in 
order to obtain the desired high-energy photons, the frequency range for the seed photons 
will change with the evolution of the electron energies, since there is a link (the limits of 
Eq. (|10|l induced by the function P(l/47^, 1, x)) between the incoming photon energy, 
the outgoing photon energy and the scattering electron energy. The electron Lorentz 
factor integration limits in Eq. 112|) have been adopted following Georganopoulos et al. 

4. Application of the model 

4.1. Parameter choice for LS 5039 

Although this model can be applied to other cases, we will test its prediction capabilities 
with LS 5039. The optical counterpart of LS 5039 is a bright (V - 11) star of ON6.5 V((f)) 
spectral type (McSwain et al. 12004(1 . at an estimated distance of 2.9 kpc (Ribo et al. 
I2()02|l . The bolometric luminosity of the companion star has been taken to be Lgtai ~ 
1.2 X 10^^ erg s~^ A new value of the orbital period, P = 4.4267 ± 0.0005 days, the 
eccentricity, e = 0.48 ± 0.06, and the semi-major axis of the orbit, a — 2.6 x 10^^ cm 
(-Rorb = o, except where otherwise stated), have been determined recently by McSwain 
et al. H2004|l . We adopt these values in our study. At radio wavelengths (1.5-15 GHz), 
the spectral index was determined to be a ~ —0.5 (where the flux density is F^, oc v") 
by Marti et al. (|1998ll . From this spectral index, we have deduced the electron power 
law index (p) following the known relationship p = \ ~ 2a — 2. From VLBI observations 
(Paredes et al. 120001 |2002|l . the radio constraint ujct cos ^ 0.15c was found. The observed 
spectrum over 100 MeV has been obtained from the third EGRET catalogue (Hartman 
et alUnSIl- The total luminosity in the range 100-1500 MeV is 2 x 10^^ erg s~^ The 
spectral energy index, rj = 0.2 ±0.2, corresponds to an observed photon index of 2.2 ±0.2 
(photon energy distribution (jj^ oc e^^'^). All these parameters are summarized in Tabled 
McSwain & Gies l|20()2|l proposed a value of 6* = 30°. Having this, and the radio 
constraint shown above {Pcos9 ^ 0.15), we obtain a wjot/c = 0.2 (Fjct = 1.02). We stress 
that our model is weakly sensitive to changes on Vjct and 9 around the quoted values 
because of the adopted low value of the Fjct. Such a low Fjet is not in disagreement 
with the value found in similar systems (e.g. 1.04 for SS 433, Spencer 119791 Hjellming & 
Johnston ll981|l . For LS 5039, we have assumed ru = 90°. In order to fix the 7-jet radius, 
we have imposed that should be neither wider than the inner disk (a few 10^ cm) 
nor thinner than some few electron Larmor radii (~ 10^ cm). Thus, has been taken 
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Parameter Description Value 



a 



orbital semi-major axis 2.6 x 10^^ cm 

Lstar star total luminosity 1.2 x 10^^ erg 

PcosO relationship between the jet velocity and the jet-observer line of sight angle 0.15 

Regret luminosity at the EGRET band 2 x 10^^ erg s"-' 

7j spectral energy index at the band 100-1000 MeV 0.2 

a spectral index at the band 1.5-15 GHz —0.5 



Table 2. Fixed parameters used in our model for LS 5039. 



Parameter 


Description 


Adopted value 


m 


angle between the jet and the orbital plane 


90° 


Pjct 


jet Lorentz factor 


1.02 


e 


angle between the jet and the observer line of sight 


30° 




7-jet radius 


10^ cm 


Lac 


accretion disk luminosity 


10"* MqC^ yr-^ 



to be 10^ cm. Regarding zq, since in this scenario (-Roib S> zq) its size does not affect 
the results, it is taken to be zero. To determine the power of the jet, we need to fix also 
the accretion luminosity of the disk (Falcke & Biermann I199(i|l . A typical value for a 
microquasar with an 0-type stellar companion can be Lac — 10"* Mqc^ yr^^. AH these 
parameters are summarized in Table |21 

To determine Lke, we will use the jet-disk coupling hypothesis for Galactic jet sources 
of Falcke & Biermann l)199()|l : Lk ~ O.l-O.OOlLac, where Lk is the total kinetic luminosity 
of the jet. Since not all the jet power is carried just by the electron population, the former 
scaling relationship turns to Lkc oc Lac , where the prior scaling factor limits are now upper 
limits. In fact, Lkc will be fixed through comparison between the observed fluxes and the 
model. From the EGRET energy range, and the involved photon and electron energy 
relationships shown above, the initial maximum Lorentz factor of the electrons should be 
about 10'*. A more accurate value will be found when we try to reproduce the observed 
spectrum slope. Regarding B^^ we will study our model behavior along a wide range of 
magnetic values, from close to estimated radio jet upper values to 100 G. 
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4.2. Results 



First, we explore the range of validity of the free parameters of the model taking into 
account observational constraints. Once this is done, we can reproduce the observational 
data with our model. Finally, additional comments on the model are made. 

4.2.1. Exploring and 755^^ 

To determine the behavior of the model with magnetic field strength, we have run our 
model for different values of the magnetic field strength and a fixed value of ^"^^ = 
5 X 10'*. We show at the top of Fig. ^ the computed IC spectral energy distribution 
[eLc) for three representative cases (1, 10 and 100 G). The spectral slopes in the energy 
range 100-1000 MeV, where EGRET detected the source, are similar. Since our model 
generates similar spectra for a wide range of magnetic field values, we have to focus 
our attention on the distances reached by the 7-jet, which is strongly dependent on B^. 
Then, high magnetic fields (tens of G or higher) would imply a 7-jet evolving very fast, 
whereas low magnetic fields (few G or lower) would imply a 7-jet evolving smoothly (see 
Sect. 4.2.3, Figs. 131 and Moreover, regarding values well below 1 G, we notice that 
upper limits of 0.2 G have been estimated from radio observations (Paredes et al. I2002|l 
for the magnetic field in the radio jet. Therefore, B-^ of about several tenths of G could 
be taken as a lower limit. 

In order to explore the dependence of the IC spectral energy distribution on 7^'^'^, 
three 7^'*'^ (10^, 5 x 10^ and 10^) have been used, fixing B^ — 10 G. The results are shown 
at the bottom of Fig.^ For 7^'*'^ = 10^, the computed 77, about 0.7, is 2.5 a softer than 
that given by EGRET observations. Therefore, 7^^'^ > 10^ is needed. An upper limit 
cannot be clearly fixed yet, since we do not know properly the spectral energy cutoff. 

4.2.2. Reproducing the EGRET data 

The computed IC spectral energy distributions for two different values of the magnetic 
field are shown in Fig. HJ We also show the EGRET data points of 3EG J1824-1514 
(Hartman et al. H^H^. The kinetic luminosity has been adopted such that the level of 
the observed flux at 100 MeV can be reproduced. The parameter values used in both 
cases are summarized in Table 13 Since we are dealing with an eccentric system, the 
star's radiation density in the jet is affected by the orbital distance variation. We have 
explored the importance of this effect by computing the IC spectral energy distribution 
at periastron and apastron passage as well as for i?orb = In the case B~^ — 10 G, 
since SSC is strongly dominant, the effect of the orbital eccentricity on the computed IC 
luminosity is not significant. However, for a B-^ = 1 G, the EC scattering is dominant 
and the orbital distance variation between apastron and periastron passages produces 
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Fig. 1. Computed IC spectral energy distributions normalized to their value at 100 MeV 
for different values of the magnetic field strength and the maximum electron Lorentz 
factor, a) For values of the magnetic field of 1, 10 and 100 G and fixing the maximum 
Lorentz factor at 5 x 10*. b) For values of the maximum Lorentz factor of the electrons 
of 1 X 10"*, 5 X 10** and 1 x 10^ and fixing the magnetic field to 10 G. 



a change of about 50 per cent in the computed IC luminosity. At the top of Fig. El we 
have plotted our results obtained at three orbital distances: 1.35 x 10^^ cm for periastron, 
2.6 X 10^^ cm for the semi-major axis of the orbit and 3.85 x 10^^ cm for apastron. Since 
the viewing periods of EGRET are longer than the orbital period of the system, we are 
not able to relate the variability predictions of our model to the variability detected by 
EGRET. As can be seen in Fig. [21 our model reproduces properly the EGRET data 
below about 1 GeV. Beyond this energy, the computed spectrum becomes significantly 
softer, but it is not in disagreement with the upper limits given by EGRET at this energy 
range. 
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Fig. 2. Two IC spectral energy distributions computed with the present model using 
the physical parameters of Table 13 The EGRET data points are also shown. The upper 
limits on undetected EGRET points are plotted with arrows, a) A magnetic field of 1 G 
has been assumed. Also, the IC spectral energy distribution for both the apastron and 
the periastron passage are shown, b) A magnetic field of 10 G has been assumed. 

Table 3. Parameter values used in the model for Fig.|2| 



Parameter 


top 


bottom 




1 G 


10 G 




10^"^ erg/s 


3 X 10=*® erg/s 


_ max 
7e0 


10^ 


10^ 



4.2.3. Further comments 

Additional comments have to be made about the magnetic field in the 7-jet. We have 
plotted the temporal evolution of the synchrotron and companion star radiation energy 
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Fig. 3. The evolution of the synchrotron and companion star radiation energy densities 
within a given shce. The value of the magnetic field is B-y = 1 G and the maximum 
Lorentz factor of the electrons is 5 x 10"'. The time axis finishes when a given slice 
moving along the jet stops emitting due to losses through the IC effect at 100 MeV. 

densities for both a of 1 G and a of 10 G within a given slice moving along the 
jet (Figs. 01 and 0J. ft shows that a variation of one order of magnitude in B^ implies 
a variation of more than two orders of magnitude in the synchrotron radiation energy 
density. Therefore, in the context of our model, magnetic fields between 1 and 10 G 
or higher could be tightly related to the high-energy 7-ray emission in LS 5039, as the 
SSC loss process is important and even dominant. Finally, we have explored the infrared 
synchrotron emission from the 7-jet. Taking a By — 10 G, we have obtained a total 
luminosity in the range lO^^-lO^** Hz of about 10'^^ erg s~^ for LS 5039. If we compare 
this value to the luminosity of the companion star in this frequency band, we find that 
the infrared jet luminosity is five orders of magnitude lower. 

5. Discussion 

As we have shown, our model is able to reproduce the high-energy 7-ray emission of the 
EGRET source BEG J1824— 1514 by assuming its origin in the microquasar LS 5039. 
Since there are no special constraints associated with a particular source in our model, 
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Fig. 4. Same as Fig. |3 but for = 10 G. 



it could be applied to other unidentified EGRET sources that might be also associated 
with microquasars. A maximum electron Lorentz factor higher than 10^ is a good pa- 
rameter range to obtain the observed spectral slope. An upper limit on 7"g'*'' has not 
been determined yet due to the lack of knowledge about the high-energy cutoff. Further 
observations by the next, more sensitive, hard 7-ray instruments (AGILE^, GLAST"^) are 
needed to determine this cutoff. To reproduce the observed luminosities, L^o ^ 10'^^ erg/s 
is needed for a magnetic field of 1 G, and Lko ^ 3 x 10'^^ erg/s is needed for a magnetic 
field of 10 G. Both values for Lkc are well within the range of allowed values established 
in the parameter discussion. 

For LS 5039, regarding the intensity of the magnetic field, the SSC effect would be 
dominant for values of about 10 G or higher, and at least non-negligible for values of 
about 1 G. In the second case, due to the importance of the star's seed photons, variations 
in the orbital distance have a significant effect on the IC flux on the scale of the orbital 
period of the system. In the case of high Bj (above 10 G), due to the strong energy 
losses, the electrons might need to be reaccelerated after leaving the 7-jet to reach the 
observed radio jet. However, for B-y ~ 1-10 G, reacceleration might not be necessary in 
regions closer to the compact object than the radio jet. 

^ http: / / agile.mi.iasf.cnr.it 
^ http://glast.gsfc.nasa.gov 
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Our approach is an attempt to join radiation and matter in a complete physical 
way. The introduction of losses is a necessary step to obtain realistic results, and a way 
to determine at the end the importance that these losses can have. The importance of 
the synchrotron process gives strong relevance to the magnetic field. Such a question is 
relevant studying the seed photon origin in the IC interaction and the timcscalcs of the 
energy processes involved, and it is also related to variability. Additional elements in our 
scenario such as the corona, disk, environment interaction, reacceleration processes, pair 
creation-annihilation phenomena and accretion variability will be considered in future 
studies. Also, emission at other wavelengths will require a proper scenario, different from 
the one presented here. 

Although future observations have to be done to confirm the nature of LS 5039 as a 
7-ray emitter, our numerical model is able to produce significant 7-ray emission levels. 
AGILE, an Italian hard 7-ray satellite expected to be in orbit in 2005, with an angular 
resolution two times better than EGRET, will make it possible to improve the accuracy 
of the hard 7-ray source position. Also, its sensitivity, almost three times better than 
that of EGRET, could define more clearly the spectral shape and its cutoff. The next 
generation of hard 7-ray satellites represented by GLAST, which will be put in orbit 
around 2008, will improve even more the study of the variability and spectral features 
at this energy range, in order to better understand the physics underlying the very high 
energy processes. 
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